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1. Introduction 

The Laplace and Fourier transforms are powerful tools widely used in scientific fields 
such as mathematics, physics, biology and chemistry. These transforms are often applied 
to linear partial differential equations and integro-differential equations to eliminate time 
and space dependence. The analytical solutions thus obtained need to be inverted to the 
time and space domain (see [1] for a monograph on the Laplace transform). Santos [2] 
found a procedure for an analytical inversion of the Laplace transform reducing the 
inversion formula to an integration on the interval [0, oo). While literature on the 
numerical inversion of the Laplace transform is rich (see [3] for a review) the analytical 
inversion still rests mostly on the tauberian theorem. This paper focuses on first and 
second order integro-differential equations containing an arbitrary kernel. It is virtually 
impossible to list a complete bibliography on the topic and for this reason the reader 
is referred to the following exemplary papers jU EJ El [7J [8] . In these works the authors 
discuss the first order integro-differential equation of the form 



d r* 

j-F(t) - zF{t) = / K(t - t')F{t')dt' 
at J 

or, after applying the Laplace transform, the equivalent equation 



(i) 



u — z — li [U) 

where, by definition F(u) = / °° exp[— ut]F(t)dt. The above equation, Eq. (j5J), represents 
a typical form that leads the process for finding a solution of the original problem stated 
by Eq. (JTJ. Once the transformed function is obtained the inversion process is often a 
difficult task. In this paper we shall consider the inversion of a Laplace transform of the 
form 

F(u) = (3) 

P n (u)+eK(u) 

and analogously, of a Fourier transform of the form 

P n {u)+eK(uj) 

where P n {u) or P n (u) is a ra-degree polynomial, e is a parameter, and K{u) or K(u) is an 
arbitrary function. Without loss of generality we shall consider the case of the Laplace 
transform. With slight changes, the results may be applied to the Fourier transform. 

The main goal of this paper is to give a prescription to find an asymptotic expression 
for the function F(t) in the representation of the starting variable, typically the time 
domain for the Laplace transform, or the space domain for the Fourier transform. This 
problem is partially solved by the use of the tauberian theorem but, as it is well known, 
the conditions for correctly applying this theorem are quite strict. We shall focus on 
the case of P\{u) = u + a and i^^) = u 2 ±a 2 giving sufficient conditions on K{u) that 



3 



allow us to find an approximate expression for either the inverse Fourier or the Laplace 
transforms. 

This work is organized as follows: In Sec. |2] we consider a short review for the 
case when the function K{u) is a generic polynomial such that an analytical expression 
for the asymptotic solution F[t) is given when P n (u) is both a first and second degree 
polynomial. In Sec. |3]we adopt the multi-scale approach to find an approximate solution 
for the case when P n (u) is a second degree polynomial and K(u) is a generic function. 
In Sec. H] the problem is solved in a less generic way but it is mathematically rigorous. 
An equation for the asymptotic expression of F(t) generated by Eq. (pQ) shall be found. 
Such an equation is independent of the form of the kernel. This is why we can use the 
term universality to describe the asymptotic equation. Finally, in Sec [5] the previous 
results are applied to the case of the convoluted Lindblad equation [9j [10] whereby we 
discuss the sufficient conditions on the kernel to obtain a positive definite quantum 
density matrix. 



2. Laplace transform containing polynomials 



In this section we briefly examine the polynomial case, however before exposing the 
main idea let us clarify a key point. Intuitively, we could say that for e — > the 
inverse Laplace transform of Eq. ([3]) can be evaluated at the poles of the unperturbed 
polynomial. Moreover, we could try to apply the tauberian theorem to "guess" the 
asymptotic solution. The following example shows that, in general, the problem can be 
much more complex. To illustrate the main idea let us consider the following Laplace 
transform 

1 



F ( u ) = ; r~, 5T- ( 5 ) 

v ' u + 1 + e O 3 + u 2 ) 

Naively we could say that for e — > 0, then F(t) « exp[— £]. First, let us find the 
solution neglecting u 3 employing the tauberian theorem idea. The poles can be evaluated 
analytically and its expressions are 

-1 - y/Y^Te 



Ui = 

Consequently the solution is 

2e~i sinh 
F{t) = - 



2s 



U-2 



-i + y/r^4i 

2e 



2e 



exp 



-t 

— + t 

e 



which seems to support the idea that for e — > 0, then F(t) 
the exact poles 

m = -l, u 2 -, u 3 



+exp[-(l+e)t] «exp[-t](6) 
exp[— t}. If we evaluate 



the exact solution is 



F(t) 



— cos 



t 



+ y/e sin 
lTe 



t 



(7) 



4 



Note that the limit for e — > of solution ([7]) does not exist. The previous example 
clarifies an important point. In general, for e — > 0, it is not correct to invert the Laplace 
transform evaluating the approximate poles of the unperturbed polynomial [see Eq. (jBJ)]. 
The reason why the expansion in power of the variable u fails in the polynomial case is 
because by neglecting higher powers we are neglecting poles that are divergent for e — > 0. 

We now consider a second degree polynomial for P n (u), namely P2{u), and at the 
end of this section we shall consider a first degree polynomial, Pi (it). Without loss of 
generality, we shall focus on the case P2{u) = u 2 ± a 2 . We start by considering K(u) = 
G n (u) where G n (u) = 'Y^l = Qa>kU k is a polynomial of n degree with n > 2. If we want an 
approximate solution of e order we must evaluate all the poles in addition to the two 
given by 



ui : 2 = ±y -a 2 - eG n (V-a 2 ). 

To fix the ideas, we select P2{u) = u 2 + a 2 . Looking for a scaling such that the term u n 
is of the same order of u 2 , we perform the transformation u — > e u U so that we have 

Equating the exponents we find that v = —^Z2- Keeping only the lowest order, we 
obtain the poles for the polynomial equation 

«n^ n + Ul = 0. (9) 
The solution can be easily found as 

iA) Ohnri 

fc = 0,---n-3 (10) 



U (k) = ^-^exp 



id) 2kiri 

-—- H 

n-2 n-2 



with = arg(— a n ). Eq. (flOj) gives n — 2 solutions that, combined with the two solutions 
given by 



"1,2 



±V-a 2 -eG n {\f^?) 



complete the n solutions for the total polynomial. The next order for the n — 2 divergent 
solutions is given by: 

Ui(k) = - r ^-. (11) 

Using the residue theorem we find the inverse Laplace transform by evaluating the 
integral 

joo+7 

F(i) * f °PM A, (12) 
2m J P 2 (m) + eG n (u) 

— 100+7 

at the poles given by 

u k = - an ~^ — , k = 0--n-3, u h2 = ±\/-a 2 ~ eG n (v / -a 2 ).(13) 

e —z (n - 2)a n v 
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The approximate expression for the function F(t) is 

exp[wit] exp[w 2 *] 



n— 3 



F(t) 



P^+eG^uJ P>(u 2 ) + eG' n (u 2 ) 



k=0 



P^(u k ) + eG' n {u k ) 



where P 2 {u) and G' n {u) are the derivatives of the polynomials evaluated in u k and u\ >2 - 
Similarly, for the poles of the first degree polynomial case, P\{u) = u + a, we obtain the 
following expressions, 

Ui = —a — eG n (—a) 

and 



Uk 



i zr ex P 

n— 1 £ n— 1 



1 



+ 



2km 



n 



1 



+ 



n 



1 



for A; = 0, • ■ • n — 2. 



3. Laplace transform containing a generic function 

In this section we shall consider a Laplace transform containing a generic function K(u) 
that can be developed in the Taylor series at the unperturbed poles given by the zeros 
of the polynomial P n {u). As shown in Sec. [2] we can not neglect the higher terms of u 
powers. Nevertheless, the asymptotic behavior of K(u) gives us information about the 
"effective" power of K{u). We set the condition that K{u) grows slower than P n (u), 
more precisely 

Km |^ = 0. (15) 

It is important to stress that the limit is performed on the positive real axis and not on 
the complex plane. 

Next we choose the case of P 2 {u) = u 2 ± a 2 which results in a typical expression, 
for example the calculus of Green's function (see Eq. (J3J), Ref. [TT] for quantum 
applications). Our goal is to evaluate the integrals 

oo 

1 p Ut f IUJX J 

F (t) = J_ / e - , du, F(x) = / s — (16) 

2m J r u 2 ±a 2 + eK{u) J u 2 ± a 2 + eK{u) 2vr 

— oo 

representing the inverse of the Laplace and Fourier transforms, respectively. In general 
the straightforward expansion of the integrand function in e powers is not correct at 
t-scale or x-scale of the order of 1/e. In other words, the expansion in e powers such as 

oo oo 

r e lu)X du f e luJX K{u) du 
u 2 ± a 2 2^ ~ 6 J (u 2 ±a 2 ) 2 2^ ( ' 

— oo — oo 

leads to an unsatisfactory approximation. It is also worthy to stress that the second 
integral on the right side of Eq. (fT7|) can be as difficult to evaluate as the initial one, 
Eq. (USD- 
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In the time representation we write the equation for F(t) as 



d 2 



df 2+a+£K [dt 



d 



F(t) = 



where the term eK f descends directly per the hypothesis that K(u) can be developed 
at the unperturbed poles in the Taylor series. Since the hypothesis states that K(u) 
grows slower than u 2 , we deduce that when t — > then F(t) « t. Using this information 
we find that F(0) = and F'(0) = 1. With regard to Eq. (PD the "order" of the 
derivative given by K (jy is smaller then the second derivative so that we can consider 
the term eK (|) F(t) slow varying function of t. Applying the multiple scale 

technique [12] we write F(t) as 



F(t) = F(r ,r l ,r 2 , 



F Q (t) + eF^t) + 



d 
dt 



5 



9 



3 



<9r + £ dri + £ dr 2 



(19) 
(20) 



where, by definition r = t, T\ — et, r 2 = e 2 t, ■ ■ -. The solution of Eq. (115]) at zero order 
in e is 



F (t) = A(n, r 2 , • ■ ■) exp[zar ] + B(ri, r 2 , • ■ ■) exp[-«ar ]. 



(21) 



To determine the functions A(ri, r 2 , ■ ■ •) and B(ti, r 2 , • ■ ■) we need the equation of the 
first order in e, 



dr 2 



+ a 1 



F 1 (t) + 2 



d 2 



drodn 



F (t) + K 



_d_ 

dr 



F (t) = 0. 



(22) 



The solution for A(t\, r 2 , • ■ ■) and B(t±, r 2 , ■ ■ •) with the conditions F(0) = and F'(0) 
1 is 



F(t) 



Now let us consider the following example with K(u) = u c 



exp 


\a(l + e^)t_ 


— exp 




2ia 


l + e 


/ K(ia) . K(-ta) \ ' 
\ Aa? 1 4a 2 J 





(23) 



F(u) 



1 



u 2 + 1 + eu c 



According to Eq. (123"]) . the solution is approximated by 

exp [-ft sin ?f] sin [t (l + f cos ^)] 



F(t) 



1 + | COS 22 



(24) 



(25) 



The numerical check is shown in Fig. [T] 

To emphasize the fact that the limit (|15p is performed on the positive real axis we 
consider K(u) = exp [—bu] as kernel in the next example. In general the condition f[T5"]) 



Figure 1. The dots represent the numerical evaluation of the inverse Laplace 
transform of Eq. (1241) . The continuous line is the approximate formula, Eq. (|25p. 
The values of the parameters are a = \/2 and e = 0.1. 



is not satisfied for «eC, with C as the complex plane. Writing F(u) with the selected 
K{u) we have 

F{u) = — — r— T . (26) 

Applying Eq. ( 12"5|) . we obtain 

exp r^te sin &1 sin \t (l + ^ecosbVi , 

Fit) = — ^ ^ — 2 - ^. 27 

w 1 + § cos b v 7 



Note that Eq. ( 127]) indicates that there is a critical value, b — 7r, for the parameter 6. This 
critical value corresponds to an exponentially growing solution or to an exponentially 
damped solution according to whether b < it or b > tt, respectively. The results are 
numerically checked in Fig. [2] 



Figure 2. Left graphic: The dots represent the numerical evaluation of the inverse 
Laplace transform of Eq. (|26p . The continuous line is the approximate formula, 
Eq. (I27p . The values of the parameters arc b = ir~ 0.5 and e = 0.1. Right graphic: The 
dots represent the numerical evaluation of the inverse Laplace transform of Eq. (I2GD . 
The continuous line is the approximate formula, Eq. (|27|) . The values of the parameters 
are b = tt + 0.5 and e = 0.1. 



4. Universality of the asymptotic equation 



In this section we will find an asymptotic expression for the inverse Laplace transform 
of a function in the form ([3]). In the previous section we focused our attention on the 
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case P2(u), whereas here we will predominantly consider the case Pi(u). This case has 
not yet been treated since the multi-scale method essentially applies to second degree 
polynomials. It shall be evident that, with slight changes, the procedure can also be 
applied to the Fourier transform. We begin with the following relatively simple case 

F(») = ^ (28) 
u — z 

where z is a complex parameter and K(u) is the Laplace transform of K(t). Making 
the hypothesis that the Laplace transform is evaluated in z, that is to say K(z) exists, 
we may write the solution of Eq. (j28p as 

F(t)= [ K(t')exp[z(t-t')]dt' = exp[zt] [ K{t') exp[-zt']dt' (29) 
Jo Jo 

and for t — > oo 

F(t) = exp[zt]K(z). (30) 

Following the same idea we now consider a more complex form of F(u). As said in 
Sec. [U a Laplace transform of the form 

F(u) = (31) 

u — z + eK [u) 

has several physical applications. As we previously assumed the function K(u) has 
an inverse Laplace transform, K(t), and K(z) exists. Considering the time domain, 
Eq. (13"T|) is equivalent to the following equation 

-F{t) - zF{i) = -e / K(t - t')F{t')dt'. 
at Jq 

Let us now examine Eq. fl3TT) in more detail. Developing its denominator in s power we 
have 

f(M) ^f ( _ ir £K. (32 ) 

w u-z ^ K ' (u- z) n K ' 

n=0 y ' 

Inverting the Laplace transform we obtain the following expression 



F W= ES^ / e n K n (t')(t~trexp[z(t-t')} = 

00 / i\n an r r t 

= J^ —r^fTn- exp ^ / K n (t')ex P [-zt'}dt' (33) 

n n - C>Z Jo 



n=0 "" "~ L J ° 

where K n {t) is the nth convolution of K(t), and in the recursive form, it is 

r-t 



K n (t) = [ K(t- t')K n _ x {t')dt', K (t) = 5(t). 
Jo 
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Since we are interested in the asymptotic behavior, we take the limit for t — > oo resulting 
in 



n=0 

oo 



n\ dz 



J2 { -^e n d " 



n=0 



n\ dz r 



exp[zt] / K n (t') exp[-zt']dt' 
Jo 

exp[zt]K n (z] 



(34) 



where, by definition, K n (z) = [K(z)] n . Applying simple algebra we obtain the following 
asymptotic expression 



F(t) = exp[zt]J2 



n=0 



(-l)" c n 

n\ 



t + 



d_ 

dz 



K n (z) 



= exp[zt](pi(e , t, z) (35) 



where the function <fii(e,t,z) is defined as 



(e,t, z) 



E 

n=0 



t + 



d_ 

dz 



K n (z) 



(36) 



Note that the parameter e is not necessarily small. The only requirement is that Eq. ([36 
has to be a convergent series. The (f>i(s,t,z) function satisfies the following equation 



9 ± , 
e—(pi{£,t,z) 



' d d 2 



(e,t,z) 



dt dzdt 

Changing the variables, t = et, v = e and w = z, leads to a simplified version 



d 2 



-(f)l(v,W,T). 



(37) 



(38) 



dwdr 

We can further transform Eq. (I38p . Considering 4>i(v, w, r) as the Laplace transform of 
a function with respect to the variable v, and making the change of variables w = x + y 
and r = x — y, we may write Eq. ( 158"]) as 

Q2 Q2 



- \§i(s } x,y) 
where, by definition 



dx 2 dy 1 



§i(s,x,y) 



(39) 



7+ioo 



$>i{s,x,y) = ^- y 0i(w,x,?/)exp[sf]^. 



7— iOO 



In Eq. f p9|) we recognize the well known Klein-Gordon equation. Using the complex 
change of variables w = x + iy and r = x — iy, we can rewrite Eq. ( 1391) as 



A$i(s,x,?/) 



9x 2 dy 2 



(40) 



that is the Helmholtz equation. The above equations; Eqs. (1381) . (|39j) and (HOI) , do 
not contain the arbitrary function ^(^). This fact implies that all integro-differential 
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equations generated by equations like Eq. (QQ) are driven by the same asymptotic 
equation. In this sense we can say that the asymptotic equation for the first order 
integro-differential equations is universal. 

To find a more manageable expression of 4>i(e, t, z), first we will rewrite it as 



E 

n=0 



-e n t r 



nl 



, Id 

1+ tlTz 



E 

n=0 



'_!)« 

v } -e n t n exp 



nl 



n log ( 1 + 



ld_ 

tdz 



K n (z) 



k{z) 



00 ( 1 ) n 



n=0 



nl 



n d 










'K(z) 


ldz_ 





00 

E 

n=0 



-l) n 

, J_ e n t n 



nl 



.(41) 



Under the condition that e is small enough to ensure a fast convergence of the series, in 
such way that only the terms n t contribute, we may write the following simplified 
expression 



00 (— l) n 

( £ ,t,z)^J2—r £nt 

Z— / 7?. I 



n=0 



cxp 



-eK(z)t 



This implies 



exp 



(42) 



(43) 



and we can call this the naive solution of the problem, that is to say, the evaluation of the 
pole in Eq. f l3T|) at the unperturbed value u — z. We stress that in general the behavior 
of 0i (e, t, z) as function of time is more complex than an exponential and the steps that 
lead from Eq. (!36|) to Eq. (JJ2J) confirm this statement. As an instructive example we 
may evaluate Eq. ( 1361) in the case where K (z) = exp[—(3z]. It is straightforward to 
obtain for <fii(e,t,z) the expression 



(e,t,z) = 



■iy 



n=0 



nl 



■e n [t-(3n] n exp [-finz] 



The series is absolutely convergent if the inequality |e/3exp [—(3z + 1]|< 1 is satisfied. 

Following the lines previously expounded for the first degree polynomial we shall 
find an expression for the case of a second degree polynomial. Let us consider the 
following Laplace transform 

1 



F(u) 



(44) 



v? ± a? + eK(u) 

For the sake of simplicity we assume that a is a real parameter and we consider the plus 
sign in the denominator. We also need the following result 



c 



-1 



1 

(u 2 + a 2 ) n + l 



1 n 



(n + k)\ 

(2a) 2 " n\ ^k\{n - k)\ 



{2at) n - k x 



7T 



at — —{n — k) 



x sin 

Following the procedures set forth in this section we find the exact expression for F 



F(t) = Im 



^ (-£)" 1 A (n + k)\ 



n=0 



a(2a) 2n n! ^ - fe)! 

x f l-2ta{t-t')] n ~ k explta(t-t')}K n (t')dt 
Jo 

and consequently the asymptotic expression 

(2iae) 



F(t) « Im 



E 



q a (2a) 2ri n! ^ ~ k)\ dz n ~ k 



Im 



E(2zae) n exp[za£] (n + k)\ 
a(2a) 2n n\ ^ 



n=0 



fc=0 



fc!(n-Jfe)! 



' d ' 


n—k 




n 






~k{z)_ 




dz 






z=ia 



Defining the function 



^ (2iaey 1 A (n + A;)! 



* 9 " 


n—k 








~k(z)_ 


dz 







r!A] ^2af n n\^k\{n-k)\ 
we may write Eq. ( I47p in a more concise way 

F{t) = Im [exp[zat]0 2 (£, a)] ■ 
As before, a further approximated expression for F(t) is given by 



F(t) « Im 



Im 



r (2met) n 
exp[zatj >^ 



^-^ a (2a) 2n n! 

71=0 V ' 



exp [za 



(2zaet)' 



, 1 9 
1 + -t0~z 



Tl 

K(- + la) 



K{z) 



^ a (2a) 2n n! 

Finally, when e is small enough, we rediscover the naive solution 



F{t) « Im 



exp 



exp 



- exp 
2ia 



-za [ 1 + 
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that basically coincides with Eq. (I2"3"j) . As for (f>i(e, t, z), the function <p 2 (s, t, a) does not 
merely represent an exponential correction to the exponential unperturbed solution. As 
previously stated, it is more complex. Note that the exponential correction holds true 
only for a sufficiently small e. 

5. Physical applications: the convoluted Lindblad equation 

A wide range of applications exists for first and second order integro-differential 
equations. In this section we shall study a case that is well known in the scientific 
literature for the quantum density matrix, the celebrated Lindblad equation [9]. We 
begin with 

^p(t) = ~[H,p(t))-L D p(t) (50) 
where Lr> is a positive semidefinite Lindblad operator of the form [10] 

LdP = —[q, [q,p]] 

with q being the systems variable measured by the environment, while 7~d represents 
the time scale of the environment-induced measurement. The Hamiltonian H and the 
observable q are properties of the system of interest, and they are operators, as prescribed 
by quantum mechanics. To reduce the number of parameters, we will set h = 1. For 
systems driven by a time independent Hamiltonian, Eq. (1501) can be promptly solved 

n 

p(t) = exp[Lt]p(0) = ^2 c k exp[\ k t]p k (51) 

fc=i 

where L = — i[H, •] — Lo is the total liouvillian, pk are its eigenvectors and its 
eigenvalues. If the matrix is an infinite dimensional matrix the sum can be extended to 
the infinite. A natural generalization of Eq. ( 1501 is 

lp(t) = f\(t-t>)Lp(t')dt>. (52) 

Barnett and Stenholm [13] studied the above equation considering the system as a 
harmonic oscillator embedded in a reservoir. They showed that, even utilizing a simple 
exponential kernel, the density matrix is positive definite only for a short time. Wilkie 
previously [JJ] showed that adding to Eq. flo^j) an inhomogeneous term with well defined 
characteristic, the positivity properties of p{t) are preserved. Recently Bologna et al. [T5] 
demonstrated that it is always possible to build a positive definite density matrix using 
the discrete version of the Lindblad equation as a starting point. The passage from 
natural time n to continuous time t is obtained performing the subordination of the 
density matrix derived from the discrete equation. 

Starting from Eq. ( |52|) it is still an open problem to find conditions on the kernel $(t) 
that generate a positive definite matrix. To achieve this, we rewrite the kernel operator 
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making a dimensionless parameter e explicit, so that $(£ — t')L — > e$(t — t')L. Then 
we evaluate the Laplace transform of Eq. (1521 

= %(()) (53) 

where $(«) is the Laplace transform of $(£). When the inverse of the operator u— e&(u)L 
acts on the density matrix it can be written in terms of the eigenvalues of L, namely 

— I — p(0) = y — * — Pk . 

u-s$(u)L ^u- e$(u)\ k 

Each addend of the sum is a Laplace transform of the form of Eq. (13T1) with z = and 
the developing parameter parameter = e\ k . Note that e k is a linear function of the 
eigenvalues. In principle, using the results of Sec. HJ and when e is small enough, we 
can write the inverse Laplace transform as 

n 

P(t) = y^c k exp[e\ k $(0)t]p k . (54) 

k=l 



As proven in Ref. [10J, the density matrix given by Eq. ( 15 ip has the proper physical 
meaning as does the matrix given by Eq. (154"]) since the only difference between the two 
is the time scale factor e$(0). But, the application of the method developed in Sec. H] 
has to be carefully considered. As stressed several times, e has to be small enough, one 
must count the factors in front of it, in order to apply the naive solution. More precisely 
the following sufficient conditions apply: 

i) |$(0)|< oo. 

ii) $(0) > 0. 

iii) e |A fc |< 1 for k = 1, 2, • • •. 

The reasons for the above conditions rest on the following arguments: Condition (i) is 
due to the requirement that in general &(z) has to exist and consequently $(0) has 
to exist (see discussion in Sec. H]). Condition (ii) is required to preserve the sign at 
the exponent in the exponential function [compare Eq. ( l5Ti) with Eq. ( 1541) ]. Condition 
(iii) ensures a small developing parameter. To insure a faster numerical convergence, 
condition (iii) may be substituted with the more convenient one 

iv) e$(0) |A fe |< 1 for k = 1,2,---. 

Considering the case studied in Ref. [13] and applying the criteria of condition (iii), we 
should select a value of e sufficiently small such that e \\ k \<^ 1 Vk. The results of Ref. [T3] 
show that the eigenvalues are linearly divergent as a function of the index, viz. ~ k, 
so that condition (iii), e |Afc|^C 1, can never be satisfied and the naive solution is not 
applicable. The density matrix evaluated in Ref. [13] is not a positive definite matrix. 
We conclude that this is due to the divergent structure of its eigenvalues. 

To elucidate this last point, we can consider a finite-dimensional density matrix 
where the condition |Afc|< M with M as a finite positive number, can be satisfied. 
The simplest case is a 2 x 2 density matrix. Strictly speaking, to keep the analogy with 
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the infinite dimensional case studied in Ref. [T3] we should consider the convolution of 
the following equation 



d 1 
-p(t) = - z [a x , [<T x ,p(t)\] 



(55) 



dt' ■ t 

where the operator q has been identified with Pauli's matrix o x . It is straightforward 
to show that the convolution of the right side of Eq. (1551) with an exponential, always 
produces a positive definite density matrix. We shall test the results of Sec. 0] considering 
the convolution of the equation containing the full liouvillian used in Ref. [T5] , 



d 1 

— pit) = Lp(t) = -iu[a x , p(t)\ --[<t z , Wz,p{t)}\ 



(56) 



where a x and a z are Pauli's matrices. The convoluted version of Eq. (!56|) may be written 

as 



dt 



p(t)=ef §{t-t')Lp{t')dt = e [ exp[-j(t - t')]Lp(t')dt 
Jo Jo 



(57) 



where e is a dimensionless parameter. The above equation admits an analytical solution. 
We shall focus on the diagonal element pn(t) since P22W can be derived via the 
relation pn(t) + P22W = 1- Taking into account that the eigenvalues of the liouvillian 
operator L are 



Ax 



A A 2 

T 



2 (-1 - VI - r 2 u 2 ) 



2 (-1 + VI - r 2 u 2 ) 



T T 

and assuming as the initial condition pn(0) = 1, we have 



, A 4 = 0, 



e 2 . , 
+ T^- 1 + 



1 



VI - T 2 CU 2 



4r ? 



Vi - 



2, ,2 



r2 cosh 



ri cosh 

tr 2 



tTt 



2t 



2t 

+ 7r sinh 



+ 7r sinh 



2t 



2t 



+ 



(58) 



where we defined the quantities Ti and T 2 as 



Ti = Jt (-8e + 7 2 r - 8eVl - r 2 w 2 ) , r 2 = Jr (-8e + -f 2 r + 8eVl - r 2 tu 2 ^ . 

From Eq. ( 1581) we can deduce that, in general, pu(t) is neither a positive quantity Vt nor a 
quantity smaller than the unit Wt. Applying the method developed in Sec.Hl the analysis 
of the eigenvalues shows that there are two cases of interest: ru < 1 corresponding to 
real eigenvalues, and too > 1 corresponding to complex eigenvalues. To simplify the two 
cases, consider tuj 1 and tu ^> 1. In the first case, condition (iv) may be written as 

As , , 

— < 1 59 

7r 
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whereas in the second case it may be written as 
< 1. 



7 



(60) 



Fig. E] shows the plot of pn(t) for given values of the parameters e, to, r, 7. In the 
left graphic, condition ( 160]) is violated because leujj^ ps 2.1, and pn(t) does not have 
physical meaning. In the right graphic, on the contrary, condition ( 1601) is satisfied 
as 2eu /j ~ 0.14, and pn(t) is a positive function smaller then the unit \/t. 





Figure 3. Left graphic: The plot of pn(t), Eq. (|58|) . The value of the parameters 
are: e = 1.5, 7 = 8, r = 1 and w = 5.5 so that 2ew/7 ps 2.1. Right graphic: 
The continuous line represents the plot of pu(t), Eq. (f58|) . while the dashed line 
represents the plot of the approximate solution, Eq. (|54|) . The value of the parameters 
are: e = 0.1, 7 = 8, r = 1 and uj = 5.5 so that 2eui/"f « 0.14. 



6. Conclusions 

This paper introduced a procedure to obtain an analytical asymptotic expression for the 
solution of first and second order integro-differential equations containing an arbitrary 
kernel. We found an asymptotic expression for the correspondent inverse Laplace 
and Fourier transforms containing an arbitrary function K(u). It was shown that a 
first order integro-differential equation is asymptotically driven by an equation that is 
independent from the specific form of the kernel of the integro-differential equation. 
A general expression for the desired asymptotic solution was also given. This result 
was applied to the convoluted version of the Lindblad equation explaining why even 
a simple kernel, such as an exponential function, does not generate a positive definite 
density matrix. Sufficient conditions on the kernel function so as to generate a positive 
definite density matrix were given at the end of Sec. El Indeed, these conditions showed 
that the structure of the eigenvalues of the liouvillian operator plays a crucial role in 
determining the positivity of the density matrix. 
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